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The direction dependence of the WIMP direct detection rate provides a powerful tool for distin- 
guishing a WIMP signal from possible backgrounds. We study the the number of events required 
to discriminate a WIMP signal from an isotropic background for a detector with 2-d read-out using 
non-parametric circular statistics. We also examine the number of events needed to i) detect a devi- 
ation from rotational symmetry, due to flattening of the Milky Way halo and ii) detect a deviation 
in the mean direction due to a tidal stream. If the senses of the recoils are measured then of order 
20-70 events (depending on the plane of the 2-d read out and the detector location) will be sufficient 
to reject isotropy of the raw recoil angles at 90% confidence. If the senses can not be measured 
these number increase by roughly two orders of magnitude (compared with an increase of one order 
of magnitude for the case of full 3-d read-out). The distributions of the reduced angles, with the 
(time dependent) direction of solar motion subtracted, are far more anisotropic, however, and if the 
isotropy tests are applied to these angles then the numbers of events required are similar to the case 
of 3-d read-out. A deviation from rotational symmetry will only be detectable if the Milky Way halo 
is significantly flattened. The deviation in the mean direction due to a tidal stream is potentially 
detectable, however, depending on the density and direction of the stream. The meridian plane 
(which contains the Earth's spin axis) is, for all detector locations, the optimum read-out plane for 
rejecting isotropy. However read-out in this plane can not be used for detecting flattening of the 
Milky Way halo or a stream with direction perpendicular to the galactic plane. In these cases the 
optimum read-out plane depends on the detector location. 



PACS numbers: 95.35.-|-d 

I. INTRODUCTION 

Weakly Interacting Massive Particle (WIMP) di- 
rect detection experiments aim to directly detect non- 
baryonic dark matter via the elastic scattering of WIMPs 
on detector nuclei, and are presently reaching the sen- 
sitivity required to detect neutralinos (the lightest su- 
persymmetric particle and an excellent WIMP candi- 
date). The direction dependence of the event rate due 
to the Earth's motion T] provides a powerful WIMP 
'smoking gun'; a directional detector needs only of order 
ten events to differentiate a WIMP signal from isotropic 
backgrounds Q, |M 13- Designing a detector capable 
of measuring the directions of sub- 100 keV nuclear re- 
coils is a considerable challenge, however. Low pressure 
gas time projection chambers (TPCs), such as DRIFT 
(Directional Recoil Identification From Tracks) |^ 
and NEWAGE |2|], seem to offer the best prospects for a 
workable detector. 

In paper I we studied the number of events required 
to reject isotropy (and hence detect a WIMP signal) and 
reject rotational symmetry (and detect flattening of the 
Milky Way halo) for a range of observationally motivated 
halo models, taking into account the detector response. 
We also calculated the number of events required to de- 
tect a deviation in the mean direction from the direction 
of solar motion due to a tidal stream. We found that if 
the senses (i.e the signs) of the recoil vectors are known 
then of order ten events will be sufficient to distinguish 
a WIMP signal from an isotropic background for all of 
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the halo models considered, with the uncertainties in re- 
constructing the recoil direction only mildly increasing 
the required number of events. If the senses of the re- 
coils are not known then the number of events required 
is an order of magnitude larger, with a large variation 
between halo models, and the recoil resolution is now an 
important factor. The rotational symmetry test required 
of order a thousand events to distinguish between spher- 
ical and significantly triaxial halos, however a deviation 
of the peak recoil direction from the direction of the solar 
motion due to a tidal stream could be detected with of 
order a hundred events, regardless of whether the sense 
of the recoils is known. While technologies for 3-d TPC 
readouts with sufficient resolution to reconstruct sub-100 
keV recoils in 3-d exist 0, , the cost and technological 
challenge of scaling these up to large, low background 
WIMP detectors is considerable. Therefore in this paper 
we repeat our analysis for a detector with less complex 
2-d read-out to assess the effects this would have on the 
detection potential. Our goals are to assess the capabili- 
ties of the next generation of detectors and present anal- 
ysis techniques which can be applied to real data (taking 
into account experimental practicalities/limitations and 
the uncertainty in the underlying WIMP distribution). 

In Sec. ^ we briefly review our calculation of the nu- 
clear recoil spectrum, including the modelling of the 
Milky Way halo. In Sec. IIIII we then apply an ar- 
ray of statistical tests aimed at probing the isotropy 
(Sec. IIII^ . rotational symmetry (Sec. IIIIBl and mean 
direction (Sec. lIII 0|) of a putative WIMP directional sig- 
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nal, before concluding with discussion of our results in 
Sec. IIVI In Appendix A we outline the circular statistics 
used. 



II. CALCULATING THE NUCLEAR RECOIL 
SPECTRUM 

The nuclear recoil spectrum depends sensitively on the 
(unknown) local WIMP velocity distribution. Observa- 
tions and simulations of halos indicate that it is likely 
that the Milky Way (MW) halo is triaxial, anisotropic 
and contains substructure (see Paper I Q for discussion 
and references) and these properties can lead to interest- 
ing features in the recoil distribution spectrum 0, 0, Q ■ 
A generic feature of triaxial halo models is a flatten- 
ing of the recoil distribution towards the Galactic plane, 
so that the recoil distribution is not symmetric about 
the direction of motion of the Sun, (^0,60). WIMPs 
from a tidal stream, with velocity dispersion small com- 
pared with its bulk velocity, produce a recoil distribution 
peaked in the hemisphere whose pole points in the direc- 
tion of the stream velocity. The net (stream plus smooth 
background WIMP distribution) peak direction depends 
on the direction of the stream and the fraction of the 
local density that it contributes. 

We consider three fiducial halo models (selected from 
the twelve considered in paper I) with properties at the 
extreme/optimistic end of the range of expected proper- 
ties. Model A (1 in paper I) is the standard halo model, 
which has a Maxwellian local velocity distribution with 
velocity dispersion equal to 270kms~^. Model B (3) is 
the logarithmic ellipsoidal model which has a mul- 
tivariate gaussian velocity distribution, with shape pa- 
rameters p = 0.9, q = 0.8 (corresponding to a density 
distribution with axis ratios 1 : 0.78 : 0.48) and velocity 
anisotropy (3 = 0.4. Model C (12) is the standard halo 
model plus a tidal stream with bulk velocity, in Galactic 
co-ordinates, (—65.0, 135.0 — 250.0) km s~^ and velocity 
dispersion 30kms"^ comprising 25% of the local den- 
sity ;i,,12,]. 

We calculate the recoil distribution for each halo model 
via Monte Carlo simulation, as described in section IIB 
of Paper I Q ^ . In a realistic directional detector it will 
not be possible to measure the nuclear recoil direction 
with infinite precision due to multiple scattering and dif- 
fusion. In paper I Q we took these effects into account 
assuming a TPC detector filled with 0.05 bar CS2, a 10cm 
drift length over which a uniform drift field of IkVcm^^ 
was applied and a 200 /xm 3-d pixel readout. Extending 
these simulations to determine the angular resolution of 
a 2-d detector is however complicated; projection effects 



The recoil momentum spectrum could also be calculated analyt- 
ically using the radon transform |^ , however for our application 
Monte Carlo simulations would still be required to calculate the 
distributions of the statistics as a function of number of events. 
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FIG. 1: The distribution of the raw 2-d angles of the recoils in 
the horizon, meridian and east-west planes from top to bot- 
tom for a detector located at Boulby for the three benchmark 
halo models (from bottom to top at <^ = 0): model A stan- 
dard halo (black solid line), model B triaxial (red dashed), 
and model C standard halo plus stream (blue dotted). We 
have set cto — 10~^ pb and po ~ 0.3 GeVcm"'^ and the in- 
tegrated distributions give the total event rate per kilogram, 
per day (0.043,0.047 and 0.050 kg"May"^ respectively). 



combined with multiple scattering of the recoil will make 
the resolution a function of the primary recoil direction 
and energy. We therefore assume a 0.05bar CS2 TPC de- 
tector with a 2-d readout of perfect resolution to provide 
benchmark figures for the numbers of events required for 
the detection of a WIMP signal with a 2-d detector. The 
numbers of events we find will hence be lower limits on 
the numbers required by a real detector. We further take 
a recoil threshold of 20 keV, as below this energy the 
recoil tracks are too short for their directions to be re- 
constructed even in our model 3-d detector. 

The orientation of the 2-d readout plane is likely to be 
determined by the geometry of the lab. We therefore de- 
fine a Cartesian coordinate system fixed in the laboratory 
in which the x-axis points towards geographic north, the 
y-axis towards geographic west and the z-axis towards the 
zenith. The three simplest possible 2-d read-out planes in 
this frame are: meridian (x-z plane), horizon (x-y plane) 
and east- west (y-z plane). We measure 2-d directions by 
projecting the 3-d recoil vectors into each of these planes 
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FIG. 2: As Fig. Q for the reduced angles, with the direction 
of motion of the Sun (which is a function of time) subtracted. 



and measuring the (right-handed) azimuthal angle </) be- 
tween the projected vector and the z-axis in the merid- 
ian plane, the x-axis in the horizon plane and the z-axis 
in the east-west plane. The orientation of the read-out 
planes also depends on the latitude of the detector. We 
focus mainly on a detector at the Boulby mine, where the 
DRIFT detector is currently located, which is at a lati- 
tude of 54.5° N. Potential locations are limited by the re- 
quirement of a suitable deep underground laboratory (in 
order to shield cosmic-ray backgrounds). All of the pro- 
posed directional detector locations which we are aware 
of lie at mid- northern latitudes (e.g. Kamioka, SnoLab). 
To cover the range of possible locations we also exam- 
ine detectors located at 36.5° N (Kamioka) and 46.1° N 
(Sudbury). 

The distributions of the angles in each plane for each 
halo model are generated from the 3-d recoil distribu- 
tions in the Galactic frame by Monte Carlo simulation of 
the sidereal time dependent coordinate transform to the 
detector frame, together with the 2-d projection proce- 
dure described above. The size of the anisotropy will be 
largest if the maximum in the recoil direction distribution 
a) spends as much time as possible close to the read-out 
plane (this minimises the spread in the 2-d distribution 
caused by projection effects) and b) has minimal motion 
in (this minimises the spread due to time-averaging). 



For smooth WIMP distributions the peak recoil direction 
is the direction of Solar motion and these requirements 
are fulfilled for any plane whose normal is at 90° to the 
spin axis of the Earth, or equivalently which contains 
the Earth's spin axis (see also Ref. |0|). This is the case 
for all meridian planes and so the anisotropy should be 
largest in this plane for any location. 

In Fig. n we plot the raw 2-d angle distributions for 
the three benchmark halo models for each read-out plane 
for a detector located at Boulby. For normalisation pur- 
poses we have taken the WIMP-nucleon cross section 
and the local WIMP density to be = 10~^pb and 
Po = 0.3GeVcm~'^ respectively. As expected the peak- 
to-trough variation is largest in the meridian plane. The 
peak-to-trough variation is smallest in the horizon plane 
as, at Boulby, this plane is furthest from the Earth's spin 
axis. The standard halo and standard halo with stream 
(models A and C) have similar peak-to-trough variations, 
while the triaxial model has a smaller variation (this is 
not so obvious from the plot, as the three models have dif- 
ferent normalisations, reflecting the different event rates 
above the 20 keV threshold) . This suggests that it will be 
hardest to reject isotropy for read out in the horizon plane 
and/or the triaxial halo model. The angle distributions 
at Sudbury and Kamioka are qualitatively similar with 
the peak-to-trough variation remaining constant in the 
meridian plane and increasing (decreasing) in the hori- 
zon (east-west) plane as the detector location is moved 
South. 

A major difference from the 3-d analysis is that here 
the recoil directions cannot be transformed from the lab 
rest frame to the Galactic rest frame. In the 3-d case, the 
time dependent conversion between the lab and Galactic 
coordinate frames tends to wash out any anisotropics in 
lab backgrounds so that they are isotropic in the Galactic 
rest frame. The modulation of the mean recoil direction 
with sidereal time (e.g. Ref. Ts'l) potentially provides a 
means of checking the Galactic nature of an anisotropic 
2-d signal observed in the lab frame. However, determin- 
ing the mean direction as a function of time necessarily 
requires large quantities of data. Instead, we use the di- 
rection of solar motion projected into each plane, /i0(i), 
which is sidereal time-dependent in the lab frame, to cal- 
culate the reduced angle, 4>red, of each recoil: 



od = - Peit) , 



(1) 



where t is the sidereal time at which the recoil occurred. 
The time-dependent nature of the reduced angle trans- 
form means that any anisotropy in lab backgrounds will 
be washed out to give isotropic reduced angle distribu- 
tions. 

The reduced angle distributions for the three bench- 
mark halo models for a detector located at Boulby are 
shown in Fig. [3 and it can be seen that these peak at 
~ 0° and have a higher degree of anisotropy compared 
with the raw distributions in Fig. ^ The reduced angle 
distributions for model C (with the tidal stream) shows 
significant excesses with respect to the distributions for 
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FIG. 3; The raw angle distributions for axial data, where the 
sense of the recoil can not be measured (halo models, read-out 
planes and detector location as in Fig. 



FIG. 4: The reduced angle distributions for axial data, where 
the sense of the recoil can not be measured (halo models, 
read-out planes and detector location as in Fig.Ql. 



the two smooth halo models (at ^lod ^ 180° — 300°, 
~ 0° ± 1° and ~ 0° - 100° in the horizon, meridian and 
east- west planes respectively). The mean reduced an- 
gles are (^^d = 348°, < 0.1° and 4.2° respectively. The 
deviation of the mean reduced angle from zero remains 
zero in the meridian plane and increases (decreases) in 
the horizon (east-west) plane as the detector is moved 
South. The reduced angle distribution for halo model B 
has a small asymmetry in the horizon plane at Boulby 
(the rate at </> is smaller (larger) than that at 27r — for 
(j) < (>)7r/2). There is a larger deviation from symmetry 
in the horizon and east-west planes at Sudbury and in 
the horizon plane at Kamioka. The reduced angle distri- 
butions in the meridian plane are always symmetric, as 
this plane is perpendicular to the flattening of the halo. 

It is possible that the absolute signs of the recoil vec- 
tors (i.e. their 'sense' -I- a; or —x) may not be measurable. 
The 2-d raw and reduced angle distribution at Boulby 
in this case are plotted in Figs. 01 and 0] The peak-to- 
trough variations of both the raw and reduced angles are 
significantly smaller than the corresponding vector angle 
distributions. The anisotropy of the raw axial angle dis- 
tributions are very small and in particular the direction 
of the stream in model C is such as to produce an almost 
flat raw angle distribution in the horizon plane. This is 



also the case at Kamioka in the east-west plane. 

III. APPLYING STATISTICAL TESTS 

As we discussed more extensively in Paper I Q a 
WIMP search strategy with a directional detector can 
be divided into three regimes i) search (aiming to de- 
tect a non-zero signal), ii) confirmation (does the sig- 
nal have the anisotropy expected for a Galactic signal?) 
iii) exploitation (extraction of information about the lo- 
cal WIMP velocity distribution, for instance flattening of 
the halo or the presence of a tidal stream). We therefore 
consider three simple hypotheses to test 

1. Is the recoil direction distribution miiform? 

2. Is the recoil distribution rotationally symmetric 
about the direction of solar motion? 

3. Does the mean direction deviate from the direction 
of solar motion? 

aimed at detecting a WIMP signal, flattening of the MW 
halo and the presence of a stream of WIMPs respectively. 

The recoil directions projected into the read-out plane 
constitute 2-d vectors, or, if the senses are not known. 
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undirected lines or axes, and so can equivalently be repre- 
sented as points on a circle, parameterised by their angle 
relative to some fixed point/direction. This allows us to 
use statistical inference methods developed for the anal- 
ysis of circularly distributed data (for a review of this ex- 
tensive field see the standard texts such as Refs. [l^LITfij V 
We investigate a variety of non-parametric statistics de- 
signed to test the isotropy (Sees. HTl A|) . rotational sym- 
metry (Sec. lIIIBjl and mean direction fSec. lIII Cjl for the 
benchmark halo models discussed in Sec. El 

For each halo model and readout plane we calculate 
the probability distribution of each test statistic, for a 
given number of events n, by Monte Carlo generating 10^ 
experiments each observing n recoils drawn from our cal- 
culated (p distributions. We then compare this with the 
null distribution of the statistic, under the assumption 
of isotropy/rotational symmetry/zero mean reduced an- 
gle calculated using analytic expressions where available 
and otherwise via Monte Carlo simulation. Specifically 
we calculate the rejection and acceptance factors, R and 
A, for each value of the statistic. The rejection factor is 
the probability of measuring a smaller absolute value of 
the statistic if the null hypothesis is true or equivalently 
the confidence with which the null hypothesis can be re- 
jected given that measured value of the statistic. The 
acceptance is the probability of measuring a larger abso- 
lute value of the test statistic if the alternative hypothe- 
sis is true or equivalently the fraction of experiments in 
which the alternative hypothesis is true that measure a 
larger absolute value of the test statistic and hence reject 
the null hypothesis at confidence level R. Clearly a high 
rejection factor is required to reject the null hypothesis. 
We also require a high acceptance, otherwise any one 
experiment might not be able to reject the null hypoth- 
esis at the given rejection factor or the null hypothesis 
might be erroneously rejected. We therefore find, using 
a search by bi-section, the number of events required for 
Ac = Rc = 0.9 and 0.95. For further details on this 
procedure see Appendix C of Paper I Q . 

With 2-d axial data the conventional (but not rigor- 
ously justified) procedure is to double the angles of the 2- 
d vectors and then reduce them modulo 360° (see Sec. lA II 
and Ref. 16]) before applying the statistical tests for cir- 
cular data. This procedure can be used in the case of the 
isotropy tests but not the symmetry and mean direction 
tests, as the effect of this transformation on the spread 
of the reduced angles is not straight forward. 

Throughout we assume zero background. This is a 
reasonable expectation for the next generation of ex- 
periments made from low activity materials with effi- 
cient gamma rejection and shielding, located deep un- 
derground As with the statistical tests of paper 
I, non-zero backgrounds can be incorporated into the 2- 
d tests as described in Section III of paper I. This re- 
quires a known background rate, so given that the next 
generation of directional detectors expect to have essen- 
tially zero background, we assume zero background to 
provide benchmark figures. We should also emphasise 



Halo 


(i?c,Ac 


= (0.90,0.90) 


(i?c,Ac 


= (0.95,0.95) 


Model 


Rayleigh 


Kuiper 


Watson 


Rayleigh 


Kuiper 


Watson 




Horizon plane 


A 


62 


70 


63 


92 


103 


92 


B 


70 


78 


71 


103 


115 


104 


C 


61 


70 


62 


91 


101 


91 




Meridian plane 


A 


18 


20 


17 


26 


29 


25 


B 


21 


23 


20 


30 


33 


29 


C 


16 


18 


16 


24 


26 


23 




East-west plane 


A 


28 


31 


28 


41 


46 


41 


B 


32 


36 


32 


47 


53 


47 


C 


28 


31 


28 


41 


46 


41 



TABLE I: Number of recoil events required to reject isotropy 
of recoil directions, A'iso, at 90(95)% confidence in 90(95)% 
of experiments for the Rayleigh, Kuiper and Watson statistics 
(Appendix A2) for fiducial halo models A (standard), B (tri- 
axial) and C (standard plus stream) for each read-out plane 
for a detector located at Boulby. 

that these non-parametric tests, unlike likelihood analy- 
ses, do not make any assumptions about the form of the 
recoil spectrum and can hence be applied to real data 
(where the underlying local WIMP velocity distribution 
is unknown). 

A. Tests of isotropy 

We initially focus on a detector located at Boulby. 
For the Rayleigh, Kuiper and Watson statistics described 
in Appendix lA 21 we determine the minimum number of 
events required to reject isotropy of recoil directions at 
90 (95)% confidence in 90 (95)% of experiments (i.e. for 
rejection and acceptance probabilities of Rc ~ Ac — 0.9 
and 0.95), A^iso, for all three benchmark halo models. 
The results for the raw angle distributions are tabulated 
in Tables HI and ITll for vector and axial data respectively. 

For vector data, in the horizon plane roughly 60 (90) 
events are required to reject isotropy at 90 (95)% confi- 
dence in 90 (95)% of experiments. The numbers of events 
required in the east- west and meridian planes are smaller 
by factors of roughly 2 and 3 respectively. For each read- 
out plane the Rayleigh and Watson tests are equally pow- 
erful, with the Kuiper test requiring roughly 10% more 
events. The triaxial model needs ^ 10% more events 
than the standard halo model while the tidal stream in 
model C does not effect the number of events required. 
These quantitative trends match the expectations from 
examining the angle distributions in Sec. II. We note 
that these number are a factor of 2 — 7 (depending on 
the read-out plane) larger than for a 3-d detector with 
perfect resolution 

The same trends (between halo models, read-out 



Halo 


(7?c, Ac) = (0.90,0.90) 


(Rc,Ac) = (0.95,0.95) 


Model 


Rayleigh 


Kuiper 


Watson 


Rayleigh 


Kuiper 


Watson 




Horizon plane 


A 


8 300 


9 300 


8 400 


12 000 


14 000 


12 000 


B 


10 000 


11 000 


10 000 


15 000 


16 000 


15 000 


C 


> 20000 


> 20000 


> 20000 


> 20000 


> 20000 


> 20000 




Meridian plane 


A 


1 100 


1 200 


1 100 


1 600 


1 800 


1 600 


B 


1 200 


1 300 


1 200 


1 800 


2 000 


1 800 


C 


1 700 


1 900 


1 800 


2 600 


2 800 


2 600 




East-west plane 


A 


2 100 


2 800 


2 200 


3 100 


3 500 


3 200 


B 


2 400 


2 800 


2 400 


3 600 


4 000 


3 600 


C 


5 500 


5 900 


5 400 


8 200 


8 700 


8 000 



TABLE H: Same as Table I for axial data (numbers quoted 
to 2 significant figures). 
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= (0.90, 0.90) 




= (0.95,0.95) 


Model 


Rayleigh 


Kuiper 


Watson 


Rayleigh 


Kuiper 


Watson 




Horizon plane 


A 


21 


33 


21 


30 


33 


31 


B 


24 


26 


24 


34 


38 


25 


C 


26 


29 


27 


38 


42 


39 




Meridian plane 


A 


12 


14 


13 


18 


20 


18 


B 


14 


16 


14 


20 


23 


21 


C 


13 


14 


13 


18 


20 


19 




East-west plane 


A 


17 


18 


17 


24 


27 


25 


B 


19 


21 


19 


28 


31 


28 


C 


19 


21 


19 


28 


31 


28 



TABLE HL Number of events required to reject isotropy of 
reduced angles. 



Halo 


(i?c,Ac 


= (0.90, 0.90) 


[Re, A, 


= (0.95,0.95) 


Model 


Rayleigh 


Kuiper 


Watson 


Rayleigh 


Kuiper 


Watson 




Horizon plane 


A 


310 


320 


310 


450 


460 


460 


B 


360 


400 


360 


520 


590 


530 


C 


670 


740 


670 


990 


1100 


990 




Meridian plane 


A 


160 


180 


170 


240 


270 


240 


B 


200 


220 


200 


290 


320 


290 


C 


280 


310 


280 


410 


460 


410 




East-west plane 


A 


220 


240 


220 


320 


360 


330 


B 


260 


300 


260 


380 


420 


380 


C 


310 


350 


320 


460 


510 


470 



TABLE IV: Same as Table HI for axial data (numbers quoted 
to 2 significant figures). 



planes and statistics) are seen for the case of axial data, 
however the number of events required are increased by 
roughly 2 orders of magnitude. This is significantly worse 
than the case of 3-d data, where we found an increase of a 
single order of magnitude 4] . The increase is even larger 
for model C in the horizon plane (as noted in Sec. II the 
direction of the stream in this case is such as to produce 
an almost flat axial angle distribution). 

We also apply the isotropy tests to the reduced angle 
distributions. The resulting numbers of events required 
to reject isotropy are given in Tables Hill and llVl for vector 
and axial data respectively. The variations between halo 
models and statistics are broadly the same as when the 
tests are applied to the raw angle distributions, however 
the numbers of events required are significantly smaller: 

30% in the meridian and east-west planes and a factor 
of ^ 3 in the horizon plane. The smaller numbers of 
events reflect the greater anisotropy of the reduced angle 
distributions. 
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Model 



A 


21 


20 


17 


30 


28 


25 


B 


24 


22 


19 


34 


32 


28 


C 


26 


24 


20 


38 


35 


28 



[Re, Ac) = (0.90,0.90) 



Boulby Sudbury Kamioka Boulby Sudbury Kamioka 



{Rc,Ac) = (0.95,0.95) 



Horizon plane 



Meridian plane 



A 


12 


12 


12 


18 


18 


18 


B 


14 


14 


14 


20 


21 


20 


C 


13 


13 


13 


18 


18 


18 



East-west plane 



A 


17 


19 


21 


24 


28 


30 


B 


19 


21 


24 


28 


31 


34 


C 


19 


23 


26 


28 


33 


38 



TABLE V: Number of recoil events required to reject isotropy 
of reduced angles using the Rayleigh statistic for each halo 
model for detectors located at Boulby (54.5° N), Sudbury 
(46.1° N) and Kamioka (36.5° N). 



Tables V and VI compare the numbers of events re- 
quired to reject isotropy of reduced angles using the 
Rayleigh statistic for detectors located at Boulby, Sud- 
bury and Kamioka for vectorial and axial data respec- 
tively. As expected from the angle distributions, the 
number of events is smallest (and constant) in the merid- 
ian plane and decreases (increases) on moving South for 
read-out in the horizon (east-west) plane. The, frac- 
tional, variation on moving from Boulby to Kamioka is 
larger for the axial angles, ~ 35% compared with ^ 20% 
for the vectorial angles. The same qualitative trends gen- 
erally occur for the raw angles but the variations are 
larger: a factor of ~ 2 (~ 3.4) for vectorial (axial) data. 
The one exception is model C for axial data. At Sud- 
bury isotropy can be rejected with < 20 000 events in all 
read-out planes, whereas at Kamioka (Boulby) more than 
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Halo 
Model 



A 
B 
C 



A 
B 

C 



A 
B 

C 



(Re, Ac) = (0.90,0.90) 



Boulby Sudbury Kamioka Boulby Sudbury Kamioka 



{Rc,Ac) = (0.95,0.95) 



Horizon plane 



310 
360 
670 



260 
310 
340 



220 
260 
320 



450 
520 
990 



390 
460 
490 



330 
390 
470 



Meridian plane 



160 
200 
280 



160 
200 
280 



160 
200 
280 



240 
290 
410 



240 
290 
410 



240 
290 
410 



East-west plane 



220 
260 
310 



250 
300 
330 



300 
350 
610 



320 
380 
460 



370 
440 
490 



440 
520 
900 



TABLE VI; As table V for axial data. 

> 20 000 events are required for the east-west (horizon) 
plane reflecting the close to flat raw axial angle distribu- 
tions for these location/read-out plane combinations. 

Finally, we translate the numbers of events required 
to reject isotropy with a detector located at Boulby into 
the equivalent detector exposures, E, required to observe 
this number of events. If the senses of the recoil direc- 
tions are observed, isotropy of the reduced angle distri- 
bution could be rejected at 95% confidence in 95% of 
experiments for WIMP-nucleon cross-sections down to 
(To ~ 7,4 and 6 x 10~^pb for read-out in the horizon, 
meridian and east-west planes respectively, with an ex- 
posure of i? ~ 10^ kg day (i.e. a 100 kg detector oper- 
ating for a period of 2-3 years), assuming a local WIMP 
density of po ^ O.SGeVcm"^. For a detector only ca- 
pable of measuring the recoil axes the sensitivity would 
be reduced by roughly an order of magnitude. We cau- 
tion once more that these numbers are for a detector with 
perfect recoil resolution and hence provide an upper limit 
on the sensitivity of a real detector. If the isotropy tests 
were only applied to the raw angle distributions these 
numbers would be signiflcantly larger, ~ 50 — 100% for 
vector data and an order of magnitude for axial data. 



B. Test for rotational symmetry 

We now examine the number of events which would 
be required to detect the deviation from rotational sym- 
metry, A^rot, for models B (triaxial halo) and C (stan- 
dard halo) using the Wilcoxon signed-rank statistic (Ap- 
pendix lA 3|) , which is applied to the reduced angle dis- 
tributions. For both models and all detector locations 
the reduced angle distributions in the meridian plane are 
completely symmetric. At Boulby the deviations from 
symmetry for model B in the other two planes are also 
small and > 20 000 events would be required. This is 
also the case for the east-west plane at Kamioka, but 
for the other location/read-out plane configurations the 



deviation from symmetry could be detected with of or- 
der 5 000 events. The difflculty of detecting flattening 
of the halo is not surprising as even with 3-d read-out 
and extremely flattened halo models thousands of events 
were required j^J . For model C the number of events is 
strongly dependent on the read-out plane and detector 
location and reflects the size of the mean reduced angle 
in each plane. The direction of the stream is close to 
perpendicular to the direction of Solar motion, so it is 
not surprising that the plane which is best for detecting 
isotropy (meridian) is the worst for detecting the stream. 
For the 10^ kg day exposure considered in the previous 
subsection, rotational symmetry could be rejected at 90% 
confidence down to (Jq ~ 1, 10 and > 50 x 10^^ pb for 
Po ~ 0.3GeVcm~'^ for model C with a detector located 
at Boulby. This test is not applicable to axial data. 



Halo 
Model 


(Re, Ac) = (0.90,0.90) 


{Rc,Ac) = (0.95,0.95) 


Boulby 


Sudbury 


Kamioka 


Boulby 


Sudbury 


Kamioka 




Horizon plane 


B 

C 


>20 000 
490 


4 800 
1 500 


5 700 
4 500 


>20 000 
760 


7 700 
2 300 


8 900 
6 900 




Meridian plane 


B 

C 


>20 000 
>20 000 


>20 000 
>20 000 


>20 000 
>20 000 


>20 000 
>20 000 


>20 000 
>20 000 


>20 000 
>20 000 




East-west plane 


B 

C 


>20 000 
4 400 


5 800 
1 800 


>20 000 
461 


>20 000 
6 900 


9 100 
2 800 


>20 000 
740 



TABLE VII: Number of recoil events required to reject rota- 
tional symmetry of recoil directions, A^rot, at 90(95)% con- 
fidence in 90 (95)% of experiments for the Wilcoxon signed 
rank statistic (appendix A3) for halo models B (triaxial) and 
C (standard halo plus stream) and each read-out plane and 
detector location. 



C. Test for mean direction 

Finally we use the Watson mean direction test (Ap- 
pendix to find the number of events required to de- 
tect the deviation from zero of the mean reduced angles 
for halo model C. The numbers of events required are, 
for a detector located at Boulby, ~ 2.5, 10 and > 10^ 
times those required with 3-d readout and are similar 
for the horizon plane and significantly (~ 50%) smaller 
in the east-west plane than those required by the rota- 
tional symmetry test. This shows that the mean direc- 
tion test is more powerful than the rotational symmetry 
test for detecting a tidal stream. For a 10^ kg day ex- 
posure considered in the previous subsection, zero mean 
reduced angle could be rejected at 90% confidence down 
to (To ~ 1,6 and > 50 x 10"^ pb for po ~ 0.3GeVcm"^ 
with a detector located at Boulby. This test can not be 
applied to axial data. 

We caution that model C has parameter values at the 
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Halo 
Model 


(Re, Ac) = (0.90,0.90) 


{Rc,Ac) = (0.95,0.95) 


Boulby 


Sudbury 


Kamioka 


Boulby 


Sudbury 


Kamioka 




Horizon plane 


C 


470 


980 


2 500 


710 


1 500 


3 800 




Meridian plane 


c 


>20 000 


>20 000 


>20 000 


>20 000 


>20 000 


>20 000 




East-west plane 


c 


2 700 


1 300 


490 


4 200 


1 900 


750 



TABLE VHI: Number of recoil events required to detect a 
deviation of the mean direction from the direction of solar 
motion, A^dir, at 90 (95)% confidence in 90 (95)% of experi- 
ments using the Watson mean direction test f appendix lA 4ll 
for halo model C (standard halo plus tidal stream) for each 
read-out plane and detector location. 



optimistic ends of the ranges estimated in Ref. i.e. 
high density and low velocity dispersion. A lower stream 
density and/or a higher velocity dispersion would give 
a peak recoil direction closer to the mean direction of 
motion of the Sun, and make the deviation due to the 
stream harder to detect. In general the directional de- 
tectability of cold streams of WIMPs will clearly depend 
on how much the projection of their bulk velocity into 
the read-out plane deviates from the projection of the 
direction of solar motion 



IV. DISCUSSION 

We have studied the application of non-parametric 
tests, developed for the analysis of circular data [l5lll6j| . 
to the analysis of simulated data as expected from a TPC- 
based directional WIMP detector with 2-d read-out. As 
the (energy and direction dependent) resolution of a 2-d 
directional detector has not been calculated to date we 
assume perfect resolution. Our results therefore provide 
a lower limit on the number of events required with a real 
detector. 

We found that if the senses of the recoils are known 
then between 10 and 30 events, depending on the read- 
out plane, will be required to reject isotropy of the re- 
duced (with the direction of motion of the Sun sub- 
tracted) angle distribution. If the senses are not known 
then these numbers are increased by roughly an order of 
magnitude. These numbers are broadly similar to those 
for full 3-d read-out. If the isotropy tests are applied to 
the raw angle distribution, however, these numbers in- 
crease significantly; 50% in the meridian and east-west 
planes and a factor of 3 in the horizon plane for vectorial 
data. For axial data the increase is even larger, at least 
an order of magnitude. Using the reduced angle distri- 
bution also has the advantage that the time-dependence 
of the transformation means that even anisotropic lab 
backgrounds will have isotropic reduced angle distribu- 
tions. It is therefore crucial that recoil events in a detec- 



tor with 2-d read-out are accurately time stamped and 
the reduced angles calculated and analysed. The num- 
ber of events required is always smallest in the meridian 
plane, which contains the Earth's spin axis, as for this 
plane the 2-d projection effects which reduce the size of 
the isotropy are minimised. 

After rejecting isotropy the next step would be to study 
the direction dependence of the signal and attempt to de- 
rive information about the dark matter distribution. If a 
significant fraction of the local dark matter distribution 
is in the form of a cold flow/tidal stream then the peak 
recoil distribution will deviate from the direction of solar 
motion, or in other words the mean reduced (with the di- 
rection of solar motion subtracted) angle will differ from 
zero. The size of the mean reduced angle (and hence 
the detectability of the stream) depends on the direction 
and density of the stream. As an example we consider 
a stream with bulk velocity, in Galactic co-ordinates, 
(-65.0, 135.0 - 250.0) kms^^ comprising 25% of the lo- 
cal density ^E^. The number of events required de- 
pends on the size of the mean reduced angle and hence 
the detector location and read out plane. In the meridian 
plane (which was best for rejecting isotropy) the mean re- 
duced angle is essentially zero, and the stream can not be 
detected, for all detector locations. For the other read- 
out planes the number of events ranges between 500 and 
5 000. In the horizon (east- west) plane the deviation of 
the mean reduced angle from zero decreases (increases) 
as the detector location is moved South from Boulby to 
Kamioka and the number of events required hence in- 
creases (decreases). The number of events required could 
in theory be reduced by optimising the choice of read- 
out plane, however in reality this is not feasible due to 
technical limitations and our lack of knowledge of the 
underlying WIMP distribution. 

It is also potentially interesting to look for deviations 
from rotational symmetry due to either flattening of the 
Milky Way halo or the presence of a tidal stream. Un- 
fortunately only a very significant fiattening of the Milky 
Way halo would be detectable. A tidal stream could be 
detected in this way, however the number of events re- 
quired is larger than for the mean direction test. 

In a realistic directional detector it will not be possi- 
ble to measure the nuclear recoil direction with infinite 
precision due to multiple scattering and diffusion. As 
discussed in Sec.mthe angular resolution of a 2-d detec- 
tor, which will be a function of recoil energy and primary 
direction due to projection effects, has not yet been cal- 
culated. We have therefore assumed a detector with per- 
fect resolution throughout. A rough estimate of the likely 
degradation in performance due to finite resolution can 
be obtained by examining the fraction of recoils retain- 
ing a sufficiently large (> 80%) fraction of their length 
after 2-d projection so that their direction can be recon- 
structed. For the standard halo model 57%, 63% and 
60% of the recoils, in the horizon, meridian and east- 
west planes respectively, retain > 80% of their length. 
Therefore the numbers of events required for a realistic 
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2-d detector with finite direction resolution are likely to 
be at least a factor of order 2 larger than the numbers 
we obtain. 

In summary, we have found that if the sense of the 
recoils can be measured, then the potential for detect- 
ing a WIMP signal (via its anisotropy) with a detector 
with 2-d read-out is similar to that for a detector with 
full 3-d read-out, provided the reduced angle distribution 
is utilised. If the senses of the recoils can not be mea- 
sured then the number of events required to detect the 
anisotropy of a WIMP signal is increased by an order of 
magnitude, which again is similar to the case of 3-d read 
out. We should caution, however, that this comparison 
assumes a detector with perfect resolution. The degrada- 
tion in performance due to the finite resolution of a real 
detector might be more significant for 2-d read out than 
for 3-d read out. 
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APPENDIX A: STATISTICAL TESTS FOR 
CIRCULAR DATA 

1. Introduction 

We first introduce the necessary basic definitions and 
terminology relating to circular statistics. For further 
information see the standard textbooks on this area 

2-d vectors are most easily paramcterised via their 
angles 4> (relative to some arbitrary fiducial direction). 
Given a sample of n 2-d vectors 0i, if we define 



n 



sm ( 



(Al) 
(A2) 



then the resultant length of the sum of the vectors, R, is 
given by — + S'^ and the mean direction cj) can be 
calculated via 



arctan (S/C) . 



(A3) 



adding tt if C < and 27r if 5 < 0, C > 0. 

The Cartesian coordinates of the centre of mass are 
denoted by {C, S) where C = C/n and S — S/n, and the 
mean resultant length is given by ^ = {C^ + S^Y^^. 

With axial data (i.e. unsigned lines) the standard 
procedure 0, 0| is to double the axial angles, reduce 
them modulo 360° and analyse the resulting vector data. 
There is no rigorous justification for his procedure and it 
is somewhat limited in its scope 16] . 



2. Tests of uniformity 

The simplest test for uniformity is the Rayleigh test, 
which uses the mean resultant length, R, which should 
be zero, modulo statistical fluctuations, for angles drawn 
from a uniform distribution. The modified Rayleigh [T9| 
statistic TZ* , defined as |0, |23| 



2ni?^ 



(A4) 



has the advantage of approaching its large n asymptotic 
distribution for smaller values of n than R. Under the 
null hypothesis that the distribution from which the sam- 
ple of angles is drawn is isotropic, TZ* is asymptotically 
distributed as xi with error of order Q^. This test 
is generally powerful, but is not sensitive to anisotropic 
distributions with zero mean resultant length (such as 
antipodally symmetric distributions). 

The Kuiper test [HISlllll is a variation of the well- 
known Kolmogorov-Smirnov test which measures the 
maximum deviation between the sample cumulative dis- 
tribution function (cdf) and the cdf of the uniform dis- 
tribution. The Kuiper test has the advantage of being 
invariant under cyclic transformations and equally sensi- 
tive to deviations between the cdfs over the entire range 
of (f>. As in the case of spherical (i.e. 3-d) data the mod- 
ified Kuiper statistic is defined as 



V* ^ V (n^/^ ^0.155 + 



0.24 

rT72 



(A5) 



where V is the (unmodified) Kuiper statistic 0, 

V = V+ +V- , (A6) 

and 



P+ ^ max ( C/i ) , i = l,...,n (A7) 



V — max Ut 



i~ 1 



(A8) 



where Ui = (j)i/27r and the Ui have been ordered so that 
Uj < Uj^i. An analytic expression for the asymptotic 
distribution of V* under the null hypothesis of unifor- 
mity is not available, so we calculate the null distribution 
numerically via Monte Carlo simulation. 

Another alternative test uses Watson's statistic 
which measures the mean square deviation between the 
sample cdf and the cdf of the uniform distribution. The 
modified statistic |^ is defined as 



0.1 

n 



(A9) 



where is Watson's statistic ITa, l23 which can be writ 
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ten as 



is normally distributed |2£ 



12n 



i=l 



1=1 



+ {n + l)U- 



12 ' 



(AlO) 



with U ~ Ui/n. As with the Kuiper statistic, the 

null distribution of U*^ has to be calculated numerically. 



Tests for mean direction 



A simple test for a given mean direction can be pre- 
formed using Watson's S statistic which directly mea- 
sures the deviation of the mean of the sampled angles 
{(f>) from the hypothesised mean direction (/io) f26l|. 
The statistic is defined as 



3. Tests for rotational symmetry 

The Wilcoxon signed-rank statistic |3| can be used to 
test for symmetry about a given direction /iq. The data 
is first rotated, i.e. 0, — > 0^ = 0^ — /iq, and the (f)i ordered 



so that 



< 



The rank of each | (pi \ amongst 



\(j)n\ is calculated and the test statistic W+ is given 
by the sum of the ranks corresponding to (pi > 0. Any 
(pi — should be removed from the data and the sample 
size n reduced correspondingly, and if more than one (pi 
has the same absolute value, then they should each be 
assigned the corresponding average rank. 
For n > 16 



W+ - n{n + l)/4 + 0.5 
[7i(n-H l)(2n-H l)/24]^/^ 



(All) 



sin(0 - ^o) 



(A12) 



where 



V2nR 



- ^cos2( 



-I 1/2 



(A13) 



is the sample circular standard error. For n > 25 5 is 
normally distributed under the null hypothesis that the 
sample is drawn from a distribution with mean direction 
Mo- 
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